README file for 'Replication Data for "Colliding in the shadows of giants: 
Planetesimal collisions during the growth and migration of gas giants"'
Philip J. Carter and Sarah T. Stewart, 2020,
https://doi.org/10.7910/DVN/5YQDU9, Harvard Dataverse.


This dataset contains data used to produce the figures found in the article:
"Colliding in the shadows of giants: Planetesimal collisions during the growth 
and migration of gas giants", P. J. Carter and S. T. Stewart, 2020, PSJ.

It also contains initial condition and final output files for all of the 
simulations, and code for reading this file format.

--------------------------------------------------------------------------------

Contents of this dataset:

pkdgrav_ss.py - Python code for reading PKDGRAV (ss) files
README - This readme file

initcond/ - Contains initial conditions for each simulation. ss files are 
                    detailed below
finaloutput/ - Final output files for each simulation, an ss file and an 
                    origin_bins file (detailed below)
collisions/ - Collision data for each simulation. File format detailed below
timesnaphots/ - Output ss files for time snapshots used in static figures in the 
                    article

--------------------------------------------------------------------------------

PKDGRAV ss files:

The ss files (both initial conditions and outputs) contain properties for each 
particle at a particular time. These files can be read using the included python 
code. Contains particle mass (solar masses), particle radius (au), particle 
cartesian coordinates (au), particle cartesian velocity components 
(au/year/(2*pi)), particle ID number, and particle color/type.

The file header stores the time, the number of particles, and a file type flag.
                
--------------------------------------------------------------------------------

PKDGRAV origin_bins files:

The origin_bins files give the mass weighted provenance data for each particle 
in the corresponding ss output file. The format is as follows:

Number of particles (N)
N lines giving fraction of mass of each particle that originated in each of 30 
0.4 au wide bins
                
--------------------------------------------------------------------------------

pkdgrav_ss.py:

Python code for reading PKDGRAV ss files.
Implements the ss class for a PKDGRAV ss file and associated origin bin data. 
Contains a brief example of how to use the code.
                
--------------------------------------------------------------------------------

collisions/
*_collisions.dat:

These files contain the data for each collision, one file for each simulation.
The format is as follows:

# Header with simulation name and column labels
One row for each collision specifying: time of impact (year), impact velocity 
(km/s), mass of target (solar masses), target ID, target x position (au), 
target y position (au), target z position (au), target x velocity (km/s), 
target y velocity (km/s), target z velocity (km/s), mass of projectile 
(solar masses), projectile ID, projectile x position (au), 
projectile y position (au), projectile z position (au), projectile x 
velocity (km/s), projectile y velocity (km/s), projectile z velocity (km/s), 
number of remnant bodies, mass of largest remnant (solar masses), mass of 2nd 
largest remnant (solar masses)


--------------------------------------------------------------------------------

Multiple simulations use the same initial condition files as detailed below:

growthis1e5_100k_1_ic.ss: growthis1e5_100k_1, growthis1e5_100k_2, 
                            growthis1e5_100k_hg1, growthis1e5_100k_hg2, 
                            growthis2e4_100k_1, growthis5e4_100k_1
growthis5e4_1_ic.ss : growthis5e4_1, growthis5e4_t2_1, growthis5e4_9, 
                        nogrow5e450k_1, nogrow5e450k_hg1
growthis5e4_2_ic.ss : growthis5e4_2, growthis5e4_11
growthis5e4_3_ic.ss: growthis5e4_3, growthis5e4_12
growthis5e4_5_ic.ss: growthis5e4_5, growthis5e4_10, growthis5e4_t2_2
growthis1e5_1_ic.ss: growthis1e5_1, growthis2e4_1, growthis2e4_hg1, 
                        nogrow2e420k_1, nogrow2e420k_hg1
growthis1e5_2_ic.ss: growthis1e5_2, growthis2e4_4, growthis2e4_hg2
growthis1e5_3_ic.ss: growthis1e5_3,  growthis2e4_2, growthis2e4_hg3
growthis1e5_4_ic.ss: growthis1e5_4, growthis2e4_hg4
GT15gm5e4_100k_1_ic.ss : GT15gm5e4_100k_1, GT15gm5e4_100k_hg1
GT15gm1e5_50k_1_ic.ss: GT15gm1e5_50k_1, GT15gm1e5_50k_hg1, GT15gm5e4_50k_1, 
                        GT15gm5e4_50k_hg1, GT20gm1e5_50k_hg1, GT20gm5e4_50k_hg1
GT15gm5e4_50k_hg2_ic.ss : GT15gm5e4_50k_hg2, GT20gm5e4_50k_hg2
mig10in1e5_100k_1_ic.ss: mig10in1e5_100k_1, mig10in1e5_100k_hg1, 
                            mig10in2e5_100k_1, mig10in5e4_100k_1, 
                            mig10in5e4_100k_hg1
mig10in1e5_100k_2_ic.ss: mig10in1e5_100k_2, mig10in1e5_100k_hg2, 
                            mig10in2e5_100k_2, mig10in5e4_100k_2, 
                            mig10in5e4_100k_hg2
mig10in1e5_100k_3_ic.ss: mig10in1e5_100k_3, mig10in1e5_100k_hg3, 
                            mig10in2e5_100k_3, mig10in5e4_100k_3, 
                            mig10in5e4_100k_hg3


